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ABSTRACT 

In local helioseismology numerical simulations of wave propagation are useful to model the interaction of solar waves 
with perturbations to a background solar model. However, the solution to the equations of motions include convective 
modes that can swamp the waves we are interested in. For this reason, we choose to first stabilise the background solar 
model against convection by altering the vertical pressure gradient. Here we compare the eigenmodes of our convectively 
stabilised model with a standard solar model (Model S) and find a good agreement. 
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1. Rationale 

One of the goals of local helioseismology is to interpret ob- 
servations of the solar oscillations in the vicinity of sunspots 
and active regions. Kilogauss magnetic fields near the sur- 
face of the Sun have a strong effect on solar oscillations, 
which cannot be computed easily using perturbation the- 
ory. Numerical simulations of wave propagation may be the 
only realistic way to model these effects accurately. Two ap- 
proaches have been proposed. T he first one is to solve the 
equations of radiative MHD fe.g. lRemoel et al] [2009') in or- 
der to model the magnetic structures and the (naturally ex- 
cited) seismic waves simultaneously. The second approach is 
to compute the propagation of wav es through a prescribed 
magnetised background mo del (e.g. ICallv &: Bogdaniri997l 
ICameron et"all l2007l I2008D . While the first approach is 
more realistic, the second approach is less computer inten- 
sive. 

Here, as a first step, we consider the problem of simulat- 
ing the time evolution of small-amplitude waves through a 
quiet Sun backgro und atmosphere using the Semi-spectral 
Linear MHD code ([Cameron et al.ir2007t) . In principle, any 
re alistic model of the upper layers of th e Sun, such as Model 
S (jChristensen-Dalsgaard et aI]|1996D . is convectively un- 
stable. It is however not possible to use Model S as a 
background model for an initial-value linear-wave calcu- 
lation, since the convective modes will grow exponentially 
and quickly dominate the solution. These convective modes 
must be dealt with. 



2. The model 

We impose two main conditions for constructing a quiet- 
Sun model, which can be used in our simulations: (CI) it 
must be convectively stable and (C2) it must support os- 
cillations which have solar-like eigenfrequencies and eigen- 
functions.We choose to begin from Model S and to modify 
it so that the resulting model best satisfies these conditions. 



There are various ways to achieve convective stability. 
For example, one may adjust the adiabatic exponent of 
the model. In this study, we choose to increase slightly the 
pressure gradient such that dzP > c^dzP, where the sound 
speed, c, and the density gradient, dzP, are retained from 
Model S. This guarantees that condition (CI) is met. The 
drawback is that the p-mode kinetic energy densities are 
too small near the surface compared to those of Model S 
and that the frequencies of the p modes are slightly too 
high. We can compensate for this by changing the sound 
speed by a few per ce nt near the surface (details and jus- 
tification discussed bv lSchunker et al.j (|in preparatioij ) . We 
arrive at a particular Convectively Stabilised Model (CSM), 
which has eigenfunctions that are as close as possible to 
Model S while the eigenfrequencies remain within about 
2% of Model S values. 

The vertical velocity eigenfunctions of CSM were calcu- 
lated numerically for each individual mode as described by 
Schunke r ct al. (in preparation). The 2D simulation box is 
145 Mm in width and covers heights from z = —20 Mm 
to z = 0.7 Mm (where z = is the base of the photo- 
sphere). The model imposes damping layers above and be- 
low to minimise reflection from the boundaries. Due to the 
height of the box, we are limited to modelling modes with 
radial orders up to 4. Figure [T] shows the eigenfunctions for 
both CSM and Model S for each ridge, f through to p4, at 
various frequencies. Figure [5] shows the relative difference 
between the CSM and Model S eigenfrequencies. 

In addition to studying the eigenmodes of the model, we 
wish to check that the power spectrum of oscillations com- 
pares favourably against that of the Sun. We extend the 
simulation to three dimensions (145.77 Mm x 145.77 Mm x 
27.5 Mm) and introduce random sources 100 km below the 
surface. The sources are spatially uncorrelated and have 
an auto-correlation function exp(— f^/2r^) where t is the 
correlation time lag and r = 5.5 minutes. A n additional at- 
tenuation term taken from iGizon fc BirchI (12002) is added 
in the momentum equation to model the finite lifetimes of 
the modes. The power spectrum of 8 hours of simulated 
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Fig. 1. Plots of Vz^/p as functions of height, where Vz is a vertical velocity eigenfunction and p is density. The dashed 
curve is the convectively stabilised model (CSM) and the solid curve is Model S. Each panel shows eigenfunctions for 
one particular radial order as labelled at the indicated frequencies (black, blue, red). 



data is shown in Figure |31 The dashed curves are the eigen- 
frequencies of Model S. Waves with phase speeds higher 
than cj/fc ~ 5.86 x 10^ cm s~^ (solid line) cannot expect to 
be modelled correctly because they hit the bottom of the 
simulation box. When comparing with real solar observa- 
tions, we find the agreement encouraging: for example, the 
mode asymmetries and location of the mode ridges com- 



pare well, but the mode amplitudes and line-widths could 
still be improved. 

In short, we have constructed a background solar 
model which is convectively stable, solar-like, and can 
be used in local helioseismology studies. Opportunities 
to further improve this model remain and ISchunker et al.l 
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Fig. 2. The relative difference between CSM eigenfrequen- 
cies (u) and Model S eigenfrcqucncies (los) as a function of 
wavenumber, for modes / through to pA. 




Fig. 3. The azimuthally averaged power spectrum of 
8 hours of simulation data. The overplotted dashed curves 
are the eigenfrequencies of Model S. The solid line is the 
phase speed at the bottom of the box, above which we can- 
not expect to reproduce the solar spectrum. The lower solid 
curve excludes the gravity modes below, that were intro- 
duced by stabilising the background. 



(|in preparation! ) present a second convectively stabilised 
model with improved eigenfrequencies. 
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